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Abstract 

The Schrodinger equation has the property that when changing the 
length scale by r — > ef and the energy scale by E — > E/e 2 , the shape 
of the wavefunction remains unchanged. The same re-scaling leaves 
the intrinsic hyperpolarizability (as well as higher-order hyperpolariz- 
abilities) unchanged. As such, the intrinsic hyperpolarizability is the 
best quantity for comparing molecules since it re-normalizes for trivial 
differences that are due to molecular size and energy gap. Similarly, 
the intrinsic hyperpolarizability is invariant to changes in the number 
of electrons. In this paper, we review the concept of scale invariance 
and how it can be applied to better understand the nonlinear optical 
response, which can be used to develop new paradigms for it's opti- 
mization. 

1 Introduction 

A decade ago, the fundamental limits of the nonlinear-optical susceptibilities 
were calculated. [U [2] These calculations were based on the sum rules, which 
are equations that relate transition moments and energies to each other, to 
simplify the sum-over-states (SOS) expressions. [3] Since the polarizability 
is rigourously optimized when all of the oscillator strength is concentrated 
in one transition, one might expect the same to be true for the nonlinear 
susceptibilities. However, the sum rules show that this approach leads to 
inconsistencies unless at least two excited states are used. 

The three-level anzatz (TLA) states that only three states contribute to 
the nonlinear-optical response when the hyperpolarizability is at the funda- 



mental limit. [H El [7] The TLA, in conjunction with the sum rules, was 
used in calculating the fundamental limits. While the three-level ansatz has 
not yet been rigour ously proven, it is always observed to hold when the 
nonlinear response of a molecule is optimized. Indeed, even Monte Carlo 
methods, which randomly assign transition moments and energies while en- 
forcing the sum rules, show that the largest hyperpolarizability calculated 
using the dipole-free SOS expression always obeys the three-level ansatz. [8] 
Given that this approach makes no assumptions about the form of the wave- 
functions suggests that the three-level ansatz may be universally true. 

The fundamental limits of the nonlinear susceptibilities define absolute 
standards to which one can compare the efficiency of a molecule. [9] In Sec- 
tion [21 we show that the ratio of the nonlinear susceptibility to the funda- 
mental limit - called the intrinsic n th -ordei intrinsic hyperpolarizability, is 
a scale invariant quantity that does not change when the size of a molecule 
is changed in a way that preserves the shape of the wavefunctions. It is 
also independent of the number of electrons in the system. As such, the in- 
trinsic hyperpolarizabilities remove the effects of simple scaling, and should 
therefore be the property of focus in structure-property studies. 

In Section [3l we show how an analysis using universal diagrams of a large 
collection of molecules with the intrinsic hyperpolarizability as the depen- 
dent variable leads to a deeper understanding of the intrinsic properties of 
a material that are responsible for optimizing the nonlinear response. This 
type of analysis leads to the proposal that modulation of conjugation may 
be one method for making better molecules that exceed the long-standing 
intrinsic hyperpolarizability ceiling of 0.03. 

In Section HI we show that certain universal features are shared by all 
quantum systems whose nonlinear susceptibilities are at the fundamental 
limit. For example, extensive numerical simulations suggest that an opti- 
mized molecule must have two and only two dominant excited states. Fur- 
thermore, the ratio of energies between the two excited states, E10/E20, must 
be about 0.49; and, the normalized transition moment to the first excited 
state must be 0.709. This suggests that underlying unifying principles may 
be at work, and raises new questions about why nature appears to require 
one and only one set of properties to optimize the nonlinear response. From 
a pragmatic perspective, these universal properties provide a set of criteria 
for searching databases of molecular properties to identify new structures 
that may have an enhanced nonlinear-optical response. 
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2 Scale Invariance 



The most general form of the Hamiltonian for N particles of mass m is, 
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where fit is the momentum of particle k, e is the particle's charge, c the 
speed of light, A the vector potential, V the scalar potential; and, where r^, 
Lfc and Sk are the positions, orbital angular momenta and spins of the k th 
particle. The potentials and momenta of the particles are functions of the 
positions. 

The Hamiltonian given by Equation [T] describes any N-electron system, 
such as a molecule or charges in a multiple quantum well. For example, the 
potential, V, can contain coulomb repulsion terms, such as — e 2 /|r*i — r*2 1 , 
spin interactions of the form v{\r\ — r*2|)si-s*2> spin-orbit coupling, and exter- 
nal electric fields; while the vector potential can describe interactions with 
external electric and magnetic fields. 

The N-electron Schrodinger equation is of the form, 

H4,{r 1 ,...r N )=E4>(r 1 ,...r N ). (2) 
When Equation [2] is transformed by fk — ► efk we get, 



1 N 



e r 



eA(efi, . . . 



+ e Wefi, . . . ; si, . . . ;Li, . 



^(efi,...) = e 2 Ei/j(en,...), (3) 



where we have used pk — > Pk/t 2 since in the position representation, V| — > 
V|/e 2 . Clearly, ip(eri, . . . ) - which has the same shape as ip(fi, . . . ) aside 
from being spatially compressed by a factor 1/e - is a solution of the Schrodinger 
equation with E — ► -Ee 2 , V(r) — > V(eri, . . . )e 2 , and A (f*i, ...)—> e^4 (efi, . . . ). 
Thus, compressing the potentials spatially by a factor of 1/e, re-scaling the 
energy by e 2 , and re-scaling the vector potential by e leaves the shape of the 
wavefunctions unchanged. By re-scaling, we specifically mean the transfor- 
mation of Equation [2] into Equation [3] with the associated re-scaling of the 
energies and vector potential. 

Upon re-scaling, the position and energy product r-rE remains invariant, 



or 



f ■ rE = constant, 



(4) 
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where r is defined as follows. The dipole moment, \x is defined as, 

N 

P=-e^2fi, (5) 

i 

where f{ is the position of the i th electron, and — e the electron charge. Thus, 
the position operator is defined by, 

r f =~- (6) 

The sum rules are calculated from the fact that the commutator of the 
Hamiltonian given by Equation [1] with any component of the position oper- 
ator yields (see appendix), 

[x,[x,H]} = -. (7) 
m 

The resulting sum rules follow by applying closure (Yin \ n ) ( n \ = -0 *° Equa- 
tion [Tl yielding. [TO] 



n \ ~ / 



N5 m , p , (8) 



where x mn are the n, m elements of the position operator, the energy 
of state k, and N the number of electrons. Note that the sum rules are 
consistent with Equation [4J i.e the product of the energy and square of the 
position operator is a constant if the number of electrons is fixed. If electrons 
are added to the system in a way that does not change the energies, then 
the transition moment will be proportional to the square root of the number 
of electrons. Specifically, Equation U] can be generalize to yield, 

f-fE = kN, (9) 

where A; is a constant. 

Any component of the position operator, for example, z = z-f, obeys the 
sum rules given by Equation [8l By convention, we define the x-direction to 
be along the largest diagonal component of the hyperpolarizability tensor. 
Thus, by j3 we implicitly mean j3 xxx . Furthermore, owing to Equation [6J 
we will loosely refer to x nm as the transition moment along the x-direction 
between states n and m. 

Using the sum rules and the three-level ansatz, it can be shown that the 
polarizability, hyperpolarizability, and higher-order hyperpolarizability are 
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bounded . [Tj [2| \TT\ 1 1 2 j In particular, the fundamental limit of the off-resonant 
polarizability (also called the zeroth-order hyperpolarizability) is, 



eft \ 2 N 



a ^ = {^)Wo' (10) 

the fundamental limit of the hyperpolarizability (also called the first hyper- 
polarizability) is, 



AT3/2 
^10 



(11) 



and the fundamental limit of the second second hyperpolarizability, is, 

Similarly, the transition moment to any excited state is bounded, and for 
the first excited state is given by, 



\rr I ■ / - -MAX (ia\ 

M - V 2^ = Xo1 • (13) 
Clearly, the fundamental limit of the n th hyperpolarizability is of the form 

(n) iV( n + 2 )/ 2 
^MAX E ,( 3 „+4)/2' ( 14 ^ 
^10 

where n\2lx is the fundamental limit of xq\, V^max 1S * ne fundamental limit 
of a, V^max * s * ne fundamental limit of /?, etc. 

Recall that the n th hyperpolarizability is of the form 

\x] n+2 

where [x] represents transition moments of the form x nm and [E] represents 
an energy difference of the form E m Q = E m — Eq. Thus, re-scaling the nth 
hyperpolarizability according to Equation d] and using Equation [9] yields, 

j.(n+2)/2 

„(«) = B = — AT(«+2)/2 m 

1 [EY+ 1 [£](3n+4)/2 JV • ^ 
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The intrinsic n hyperpolarizability is given by, 

„{n) / p \ (3n+4)/2 

Vint = « ( w) ' 

where we have used Equations [14] and [TBI r/^ T is clearly unchanged under 

(n) 

simple scaling and independent of the number of electrons. Thus, Tjj^ T is a 
scale-invariant quantity. 

The intrinsic hyperpolarizabilities are quantities that remove the effects 
of scaling, and can be used as a metric for comparing a molecule's nonlinear- 
optical efficiency, independent of the number of electrons or energy gap. Sim- 
ply stated, larger molecules with more electrons will generally interact more 
strongly with light than smaller, electron poor systems. The intrinsic hyper- 
polarizabilities remove such effects, allowing one to focus on the structural 
properties that affect the response. Only then can truly new paradigms be 
developed for making large molecules with exceptionally enhanced response. 



3 Using Scale- Invar iance to Assess Molecules 

As described in the previous section, the increase of the hyperpolarizability 
of a molecule with size can be partially explained by the increase in the 
number of electrons and lowering of the energy gap (defined as the difference 
in energy between the ground and first excited state). This dependence 
we shall refer to as Simple Scaling. The intrinsic hyperpolarizability is a 
quantity that accounts for Simple Scaling, and can be used to compare 
molecules of all shapes and sizes. A plot of the intrinsic hyperpolarizability 
as a function of a parameter of interest is referred to as a Universal Plot, 
and provides a bird's-eye view of the behavior of a large and diverse set of 
molecules. 

Figure [TJ shows a Universal Plot with Xmax as the independent param- 
eter. The dashed green line, sometimes referred to as the apparent limit, 
represents an intrinsic hyperpolarizability of about 0.03. They yellow points 
show the best molecules measured prior to 2007, all of which are below the 
apparent limit. While three decades of intense research led to molecules 
with large hyperpolarizabilities, their intrinsic hyperpolarizabilities showed 
little improvement. Thus, all of the research efforts prior to about 2007 that 
were aimed at improving the nonlinear-optical response essentially lead to 
larger electron-rich molecules with smaller energy gap. 

Molecules Dl through D4 have large hyperpolarizabilities, and are used 
in making electrooptic devices. [13J However, the Universal Plot clearly shows 
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Figure 1: A Universal Plot that includes all of the best measured molecules. 
With permission from NLOsource.com 
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that the intrinsic hyperpolarizabilities are unremarkable. Indeed, molecule 
D4 has a large hyperpolarizability owing mostly to its small energy gap 
(^max is the largest of the group) and large number of electrons. If it 
were possible to adjust the structure of such molecules to attain an intrinsic 
hyperpolarizability of unity, the hyperpolarizability could be improved one 
hundred fold. 

In an effort to design new paradigms for breaking through the apparent 
limit, Zhou et al used numerical optimization techniques to determine what 
potential energy functions would lead to the largest hyperpolarizability. [5j [7] 
They found that there were many one-dimensional potential energy func- 
tions that lead to intrinsic hyperpolarizabilities of 0.709, but not one was 
found to be larger. In all cases, the three-level ansatz was obeyed. 

One class of potentials that optimize the hyperpolarizability have ran- 
dom (i.e. non-periodic) undulations that appear to serve the purpose of 
spatially separating the wavefunctions to insure that the three-level ansatz 
in obeyed. [7] As such, Zhou and coworkers proposed that modulation of con- 
jugation might be a new paradigm for making better molecules. The most 
common molecules with large hyperpolarizability where electron donors and 
acceptors separated by a conjugated bridge of the form C-C=C-C-C=C- 
C=C. The concept of modulated conjugation suggests that bridges with 
differing atoms, such as C=N-C=S-C... etc., might yield a larger response. 

Perez-Moreno and coworkers showed that molecules with a modulated 
bridge, as shown by the points labeled LCW6 and LCW7 in Figure [TJ had 
larger intrinsic hyperpolarizability than the non-modulated versions, labeled 
LCW4.p~JJ [15] While it is difficult to theoretically associate an effective 
potential energy profile of these molecules (shown in Figure [2]), they argue 
that differences in aromatic stabilization energy of the rings serve as a proxy 
for modulated conjugation. 

As the molecules in the series labeled CI to C5 are made longer (see 
Figure [2]), the intrinsic hyperpolarizability at first increases, then decreases 
as shown in Figure HJ[16j The decrease is attributed to the inability of the 
larger molecules to retain planarity. To test this hypothesis, each molecule in 
the series was placed in an amylose helix, which keeps the molecules planar. 
The measured hyperpolarizabilities are labeled Cl-III to C5-III. The longest 
molecule, C5-III, when place in the amylose helix, has the largest intrinsic 
hyperpolarizability. Its absolute hyperpolarizability is enhanced by about a 
factor of ten over the non-encapsulated form. 

Kang and Coworkers considered twisted molecules, as shown in the bot- 
tom of Figure HdZKH] The points in Figure [J labeled M2, MC2, and MC3 
clearly have the largest intrinsic hyperpolarizabilities ever measured. Thus, 
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Figure [0 
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this paradigm shows great promise in making molecules with exceptional 
hyperpolarizability if the enhancement due to the twist can be maintained 
for larger forms of such a molecule. 

4 New Paradigms for Making Molecules 

Aside from the twisted molecules, the hyperpolarizabilities of most molecules 
fall far short of the fundamental limit. Even some numerical optimization 
methods seem to be bound by the apparent limit. For example, Keinan, 
Wang and coworkers use a computer algorithm to sort through a library 
containing about a million building blocks of linear combinations of atomic 
potentials to find structures with large nonlinear-optical response. |19|. [20] 
They have found calculated off-resonant hyperpolarizabilities on the order of 
(3q = 10, OOOx 10~ 30 esu; but, the intrinsic hyperpolarizability is (f^ 1 = 0.013, 
more than a factor of 2 from the apparent limit and just over 1% of the 
fundamental limit. If the intrinsic hyperpolarizability were optimized so 
that all of the electrons were used efficiently, the hyperpolarizability would 
be about /3 = 800, 000 x 10~ 30 esu. 

In contrast, our approach is to use abstract numerical optimization, 
which apply to any quantum system, with the goal of identifying common 
features that are shared by systems whose nonlinear response is at the fun- 
damental limit. While this approach may not lead to specific molecular 
structures, it can provide general guidelines. Our numerical optimization 
studies have found certain universal properties of a quantum system when 
the intrinsic hyperpolarizability is optimized. They are as follows: 

When the intrinsic hyperpolarizability P^xx" * s optimized by varying the 
potential energy function or the geometry of point charges in 2-dimensions, 

• the ratio of energies is always about E20/E10 = 0.49; 

• the transition moment to the dominant state is given by xiq/x^q AX = 
0.79; 

• the optimized intrinsic hyperpolarizability is always 0.709 for a broad 
range of potentials; 

• exactly three states dominate the nonlinear-optical response; 

• and, the electron cloud of the three dominant states tends to be aligned 
in one dimension. 
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These observation give clues about how to design better molecules, and 
introduce new puzzles. For example, searches for better materials should 
seek molecules with approximately equally-spaced energy levels and two 
dominant excited states. One design paradigm, suggested above, is the 
synthesis of new molecules that exhibit modulation of conjugation. New 
theoretical studies, on the other hand, should focus on identifying ways of 
arranging atoms in a molecule in such as way as to provide an effective 
potential along the conjugated pathway with the kinds of undulations that 
cause the eigenfunctions to be spatially separated in a way that enforces the 
three- level ansatz. 

While these universally-observed properties might be used to design bet- 
ter molecules, the fact that universal properties exist is puzzling. Why are 
there so many different potentials that all yield the same optimized fi^Jl 
Why is the optimized value given 

by filxl = 0.709, and not unity? Fi gure 
[3] shows a two-dimensional representation of the hyperpolarizability land- 
scape, where Parameter 1 and Parameter 2 quantify the potential energy 
function used to calculate the hyperpolarizability. The same kind of land- 
scape is found for other quantities, such as the energy ratio. Why do such 
a diverse range of potentials all yield the same energy ratio and normalized 
transition moment? Answers to such questions may lead to a better under- 
standing of light-matter interaction; but, may also lead to new fundamental 
principles. 



5 Conclusion 

We have shown that the fundamental limit of the n th hyperpolarizability de- 
fines an absolute standard for comparing molecules. The ratio of a quantum 
system's n th hyperpolarizability to the fundamental limit defines the intrin- 
sic n th hyperpolarizability, which is a scale-invariant quantity that can be 
used to generate a Universal Plot of the nonlinear response of a broad range 
of molecules. Since this Universal Plot takes into account Simple Scaling, 
it can be used to better focus on the relevant properties required to de- 
velop new paradigms for understanding structural properties that enhance 
the response. 

All quantum systems with optimized intrinsic hyperpolarizabilities are 
dominated by two excited states; and, no system has been observed with an 
intrinsic hyperpolarizability near unity with more than 2 dominant excited 
states. An investigation of the potential energy functions that optimizes (3 XXX 
has led to the proposal that systems with conjugated modulation may be a 
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Figure 3: The intrinsic hyperpolarizability landscape, which abstractly rep- 
resents the fact that many optimized potential energy functions lead to the 
same intrinsic hyperpolarizability. 
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new paradigm for large-hyperpolarizability molecules. Limited experimental 
evidence supports this hypothesis. Other paradigms that appear to enhance 
the hyperpolarizability are twisted molecules; and, in the other extreme, 
using an amylose helix inclusion complex to force large chromophores to be 
planar. 

For one-dimensional quantum systems that are representable with a po- 
tential energy function, the ratio of the first to second excited state energy 
is found to be about 0.49 and the normalized transition moment to the 
first excited state is 0.79. Thus, it may be possible to identify candidate 
molecules by screening their linear optical properties for these values. More 
importantly, these results hint of a set of universal properties that are re- 
quired of optimized molecules. The fact that many different systems all 
share these universal properties is a puzzle whose solution may lead to a 
deeper understand of the nonlinear-optical response of a quantum system. 

Acknowledgements: MGK thanks the National Science Foundation 
(ECCS-0756936) and Wright Paterson Air Force Base for generously sup- 
porting this work. 



6 Appendix 



In this appendix, we show that the commutator given by Equation [7] follows 
from the general Hamiltonian given by Equation [TJ We do so by considering 
each term in the Hamiltonian. 

The potential V (r%, . . . rjv; s%, . . . sjy; L\, . . . Lj^j . i.e. the last term in 
Equation [TJ can include many types of interactions. For potentials of the 
form V (fi, . . . r/v; s±, . . . sjv), [x, V] = 0. Spin-orbit coupling contributions 
have terms of the form vS ■ L, where it is understood that v is a function of 
the electron coordinates. The commutator with x yields, 



x,vS ■ L 



vS- 



x, L 



+ v 



x, S 



■L=m, 



where f(r) is a function of position. We have used the fact that 

because S and r are operators in orthogonal sub-spaces; and that 
terms of the form 



x, S 
x, L 



[x,L z ] 



-yih. 



(18) 

= 
has 

(19) 



[x, xp y - yp x \ = - [x, yp x \ = -y [x,p x \ 

Clearly, Equation [19] gives terms that are only functions of the position. For 
terms of this type, 

[x, [x, V}} = 0. (20) 



13 



Thus, Equation 1201 generally holds. However, we note that it is possible to 
pick potential functions that will yield [x, [x, V]] 7^ 0, but none of these are 
physically meaningful. 

Next we consider the commutator of first term in Equation Q] with x. 
Expanding the square yields three terms, 



7i = [x,p 2 ] = 2ihp x , 



x,p • A + A • p 



(21) 
(22) 



and 



T 3 = [x, A 2 ] = 0. (23) 
Each term in Equation [22] is a function of the vector potential, as follows, 



x,p ■ A 



x,A 



+ [x,p\ ■ A = + ihA x . 



Using Equations [21] through [25] yields 



N 



m 



(24) 



(25) 



Equations [20] and [25] lead to Equation [7] 
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